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Abstract The equation for two reggeized gluons with a running QCD coupling constant, 
proposed earlier on the basis of the bootstrap condition, is investigated by the variational 
technique in the vacuum channel. The pomeron intercept and slope are calculated as a 
function of the infrared cutoff parameter m ("the gluon mass"). The intercept depends 

weakly on m staying in the region 0.3 0.5 for values of m in the range 0.3 1.0 GeV. 

The slope goes down from 5.0 GeV~ 2 to 0.08 GeV~ 2 in this range. The optimal value of m 

seems to lie in the region 0.6 0.8 GeV with the intercept 0.35 and slope 0.35 0.15. 

The character of the pomeron singularity: pole or cut, is also discussed. 
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1. Introduction. We have recently proposed a new method to incorporate a running 
coupling constant into the equation governing the behaviour of two reggeized gluons [ 1 ]. 
It is based on the so-called "bootstrap condition" of L.Lipatov [ 2 ], which requires that 
the solution in the gluon colour channel should be the reggeized gluon itself. The bootstrap 
condition leads to a relation between the gluon interaction U and its Regge trajectory uj: 
both have to be expressed via the same function r], so that in the momentum space 

u(q) = -mUq) J (d 2 q 1 /(2n) 2 )/r,(q 1 )r,( q2 ) (1) 

and 

U(q,quq[) = ~3(ri(qi)/v(q'i) + ^M^Mtt - q'i) + Mq)/M)M) (2) 

In (2) q is the total momentum of the two interacting gluons. In both equations q\ + q 2 = q- 
The asymptotical behaviour of 77(g) at high momenta can be found from the comparison 
of the pomeron equation at q = and the standard Gribov-Lipatov-Altarelli-Parisi evolution 
equation in the double logarithmic approximation. With a running QCD coupling constant 
g(q) we obtain at high q 

1/V(q) = 9 2 (q)/W = 2Tr/bq 2 lnq 2 /A 2 (3) 

where 

b = (11 - (2/3)iV F )/4 (4) 

and A ~ 200 Mev is the standard QCD parameter. With a fixed coupling constant (the 
BFKL pomeron [ 3 ]) l/rj(q) results proprtional to q 2 . It follows that the introduction of a 
running coupling constant changes the pomeron equation already at leading order, so that it 
cannot be studied by a perturbative approach like in [ 4 ]. 

The behaviour of rj{q) at low q 2 remains unknown and, as we want to stress, hardly 
amenable to any theoretical approach, since it belongs to the nonperturbative confinement 
region. More than that, the two-gluon equation itself looses sense at low q ~ A, where gluons 
cease to exist as individual particles. Therefore the behaviour of rj(q) at low q 2 plays the role 
of a boundary condition to be supplemented to the pomeron equation, which effectively takes 
into account the confinement effects. One can at most try to parametrize it in some simple 
manner and study the consequences. In [ 1 ] we have chosen a parametrization 

l/r)(q) = 2ir/b(q 2 + m 2 ) ln{(q 2 + m 2 )/A 2 ) (5) 

where the parameter m > A has a meaning of an effective gluon mass. 

The asymptotical behaviour of both the gluon trajectory and interaction with the choice 
(5) was studied in [ 1 ]. Taking for them simple analytic expressions which interpolate 
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between low and high values of their arguments, the pomeron intercept was also studied by 
the variational technique in [ 1 ]. 

In this note we investigate the pomeron intercept with exact values for the trajectory and 
interaction, which follow from (1) and (2) with the choice (5). We confirm that the intercept 
depends very weakly on the value of the gluon mass m, as was found in [ 1 ] for approximate u 
and U. Its absolute values turn out, however, somewhat higher: for m growing from 0.3 GeV 
to 2 GeV the intercept falls from 0.5 to 0.25. 

Since we know the pomeron equation with a running coupling constant also for q / 0, we 
are in a position to calculate the slope of the pomeron trajectory at q = 0. It happily results 

positive, its values depending very strongly on m: in the same interval m = 0.3 2. GeV 

the slope goes down from 5.0 GeV~ 2 to 0.01 GeV~ 2 . From the experimental observations 
then the optimal value of m lies in the region of the p-meson mass m = 0.6 0.8 GeV. 

We finally discuss the character of the pomeron singularity in the j-plane. As it seems, 
it depends on the chosen behaviour of 77(g) at q — > 0. For the choice (5) it turns out to be a 
moving pole, although we cannot exclude that other choices may result in a moving cut. 

2. The intercept. For the forward scattering, q = 0, the pomeron equation aquires the 
form of a Shroedinger equation in two dimensions 

Hip = Eip (6) 

where the Hamiltonian is a sum 

H = T + U + Q (7) 

of a kinetic energy 

T(q) = -2u)(q) (8) 
a local interaction U(r) with a Fourier transform 

U(q) = -6/ V (q) (9) 

and a separable interaction 

Q = (1/12) V (0)\U)(U\ (10) 

The intercept A is the ground state energy of (6) with a minus sign. 

In [ 1 ] the analytic forms were used for T(q) and U(r) which interpolate between small 
and large values of their respective arguments 

T(q) =T(0) +lnln(g 2 + m 2 ), U(r) = lnln(l/r 2 + m 2 ) (11) 

(in units 3/b, A = 1). Numerical calculations with rj(q) given by (5) lead to T(q) and 
U(r) shown in Table 1. for m = 3.0 in the regions where they noticeably differ from their 
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approximate expressions (11) (shown by dashed curves). As one observes, the largest error 
in using (11) is for U(r) at large r, where the exact potential falls exponentially, whereas an 
approximate one only as 1/r 2 . Since the potential is attractive, one might think that with 
the approximate potential the ground state level should lie deeper. However the approximate 
kinetic energy lies systematically above the exact one. This effect turns out to be dominating. 

Taking exact numerical values for T(q) and U(r) we investigated the ground state energy 
Eq by the standard variational method. As in [ 1 ] trial functions were chosen as linear com- 
binations of the two-dimensional harmonic oscillator functions for zero angular momentum 

N 

i>(r 2 ) = c k L k (ar 2 ) exp(-ar 2 /2) (12) 

with L/% the Lagerre polinomials and a and c& variational parameters. Up to 16 polinomials 
have been included in the calculations. The resulting intercept is shown in Table 2. (the 
second column) as a function of m for A = 0.2 GeV and three flavours (b = 9/4). As one 
observes, it depends on m rather weakly except for the values of m quite close to A. For 

reasonable values of m in the range 0.3 1. GeV it goes down from 0.5 to 0.3. Comparing 

these values with the ones found in [ 1 ] with (11) we find that the exact trajectory and 
interaction lead to the intercept twice as high. 

3. The slope. To determine the slope at q = we can use the perturbative approach. 
For small values of q we present the Hamiltonian in the form 

H = H + W(q) (13) 

where Hq is the Hamiltonian for q = and W is the perturbation proportional to q 2 . The 
slope a' is then calculated as the average of W in the ground state ipo at q = divided by q 2 . 

At q 7^ the pomeron equation retains its form (6). The components of the Hamiltonian 
(7) are now given by 

T = -u 1 - 0J 2 (14) 

u = -3(\fvi/mV\fm/m + y/m/mVy/m/m) ( 15 ) 

Q = 3 V (q)\l/^jm)(l/^W 2 \ (16) 

Here and in the following u± = uj(qi), rj\ = rj(qi) and so on. 

Let ft = (l/2)q + I, q 2 = (l/2)q - I. To find W{q) we develop (14)-(16) in powers of q 
up to the second order. Straightforward calculations lead to 

T = T - (1/8) [qr[qr, T \] (17) 



U=Uo + (l/2)[(ql)( 1 [(ql)( 1 ,U ]} 



(18) 



5 



Q = Qo + q%Qo ~ {(V4)g 2 Ci + (l/2)((//) 2 (C2 - Ci), Qo} (19) 
where To, Uq and Qo correspond to q = and are given by (8)-(10). We have denoted 

Ci(/) = (ln/x+l)//xln / x, ( 2 (l) = l/fx 2 \nfi, ( 3 = Ci(0) 

where n = m 2 + I 2 . Note that in the double commutators (17) and (18) one can change To 
and Uq to Tq + Uq = Hq — Qq, since the added terms commute with the respective operators. 
In this manner we finally obtain 

W = -(l/8)[qr[qr,H - Q }} + (l/2)[qlCi[qKi,H - Qo}} 

+ q\sQo ~ {(V% 2 Ci + (1/2)(?0 2 (C2 " Ci 2 ),Qo} (20) 

Taking the average of (20) in the ground state ipo we present the corresponding change 
in energy AE as a sum of three terms 

3 

AE=J2 AE(n) ( 21 ) 

n=l 

Here terms with n = 1, 2 and 3 correspond to parts of (20) proportional to the ground state 
energy Eq, containing To + Uq and containing Qo respectively. 

A£« = E ((-l/4)(qr) 2 + (ql) 2 ( 2 ) (22) 

AE^ = ((1/4) ( qr ) (T + U )(qr) - (gZ)Ci(T + E/ )(<?0Ci> (23) 

and 

Ai?® = 3 % <M(V4)(?r) 2 + ^ 2 (Cs - (l/2)Ci) - (ql) 2 (2)V)(V\^ ) (24) 

where we denoted (^ol^lV'o) = (A). 

In most of these averages one can immediately integrate over the azimuthal angle sub- 
stituting (ql) 2 — > (l/2)g 2 / 2 and (qr) 2 — > (l/2)g 2 r 2 . Two averages, however, are a bit more 
complicated. 

Ai = <(5r)ro(gr)), A 2 = <(g0& W)Ci> 

The first one is evidently a momentum space average in the state (qr)-ipo(l 2 ) = 2i(ql)^' (l 2 ), 
so that A\ = 2<7 2 (?/>ol^ 2 To|V'o)- The second average, also in the momentum space, reduces to 
that of a nonlocal potential U\\ A 2 = (l/2)q 2 {Ui) where U\ possesses a kernel 

U^hM) = -Ki{h)Uh){hh)h{h-l2) (25) 
The final formula for the slope consists of three terms in correspondence with (22)- (24) 

o! = £ a' n (26) 

n=l 
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They are 

a[ = £ <(-l/8)r 2 + (l/2)/ 2 C 1 2 ) (27) 

a' 2 = ((l/8)r 2 U - (1/2)Z 2 C 1 2 T - (l/2)Lq> + (1/2)(^|/ 2 T |^) (28) 

a' 3 = 37? <M(l/8)r 2 + Cs - (l/2)Ci - (l/2)/ 2 C 2 )V> Wo> (29) 

We have calculated a' according to these formulas taking for the ground state ipo the trial 
function (12) with the parameters which minimize the energy. The results are shown in Table 
2. (the third column). They strongly depend on m, reaching quite large values for m close 
to A and rapidly falling with its growth to values of the order 0.01 GeV~~ 2 for m = 10A. 

4. Pole or cut? We finally want to comment on the character of the pomeronic singu- 
larity in the j plane for our eqation. It is easy to show that with the choice of r)(q) according 
to (5) it is a moving pole. To prove it we have to demonstrate that the ground state wave 
function ■i/'o is normalizable. 

With r)(q) given by (5), the interaction falls exponentially for large intergluon distances 
r — ► oo, so that the equation reduces to 

(T(q) - E)^{r) = (30) 

where q 2 = —A. Evidently the asymptotics of the solution is a linear combination of zero 
order Bessel functions 

V>(r) ~ aR ( o ] (qr) + 6h[, 2) (qr) (31) 
where q is determined from the equation 

T(q) = E (32) 

However, for positive q 2 the kinetic energy T(q) is greater than a positive threshold value T(0) 
(Table. 1). Then for E > T(0) the asymptotics (31) is oscillating and ip is not normalizable. 
For E < T(0), and for E < in particular, Eq. (32) gives pairs of complex conjugate solutions 
q = qi ± q2 with q 2 > 0. Then evidently the function i/j(r) becomes exponentially falling at 
large r like exp(— q 2 r) 

To show its normilizability we have additionally to study its behaviour at r = 0, that is, 
at large momenta q. With the asymptotics of the kinetic energy clear from (11), the equation 
at q — > oo becomes 

\n\nq 2 il>{q)+bm 2 \nm 2 /q 2 \nq 2 = (l/n) j d 2 q'ij(q')/({q - q') 2 + m 2 ) ln((q - q') 2 + m 2 ) (33) 
Two alternatives are possible: either 

^{r = 0) = [ (d 2 q/2ir)ip(q) < oo (34) 
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or ip(r = 0) = oo. We want to show that the first one is realized. In fact, assume that 
ip(r = 0) = oo. Introduce a function /(£) = g 2 lng 2 ?/;(g) with £ = lnlng 2 . By assumption 

J de/(o = °° 

After calculating the asymptotical behaviour of the integral term in (33) one arrives at an 
equation for / of the form 

which can never be satisfied, since the second term on the righthand side exactly cancels the 
lefthand side term. Thus (34) is true. Then the integral term in (33) behaves like l/q 2 lng 2 
at high q, which contribution should be cancelled by the one from the separable interaction 
(the second term on the lefthand side). The first term, coming from the kinetic energy, then 
represents a subdominant correction. With (34) fulfilled, the function ip(r) is well-behaved 
for all r and is normalizable. Then solutions for Eq < should belong to a discrete spectrum. 

We have to point out, however, that this conclusion depends heavily on our choice of 
rj(q), which assumes that it is different from zero at q = 0. As discussed in the beginning, 
actually we do not know anything about the behaviour of r)(q) at q << A nor does it make 
any physical sense. More general parametrizations, like 

77(g) = (6/2vr) (g 2 + m\) ln(g 2 + m 2 ) 

with two different masses m and m\ are equally admissible. In particular we can take ra\ = 
and choose 

v (q) = (6/2vr)g 2 ln(g 2 + m 2 ) (35) 

when 77(0) = 0. This choice seems to have some advantages: now the trajectory will go 
through 1 at q = ( for the "physical" gluon) and also the solution will satisfy the condition 
ip(q = 0) = 0, which follows from the colour current conservation if one goes down to q = 0. 
However these advantages are completely spurious: with confinement one can never go as far 
as q = and should always stay at momenta q > A in coloured channels. So, in our opinion, 
the choice (35) is nothing better than our earlier (5). 

It however leads to a considerably more complicated gluonic Hamiltonian. With (35) the 
kinetic energy and interaction do not exist separately, but have to be taken together in the 
integral kernel exactly in the same manner as with the BFKL equation [ 3 ]. The infrared 
behavior of this new Hamiltonian is completely different from the one with (5) and becomes 
similar to that of the BFKL Hamiltonian. With that, the equation with a running coupling 
constant does not possess the scaling symmetry and hardly admits an analytic study. We 
cannot say anything definite about the pomeron singularity for it: it might well be a cut as 
for the BFKL equation (only a moving cut, contrary to the BFKL case). 
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Thus we conclude that in the real world with a running coupling constant the equation for 
the pomeron is sensitive to the confinement effects, that is, to the behaviour of the trajectory 
and interaction in the infrared region. The dependence on this behaviour is minimal (and 
quite weak, indeed) for the intercept, whose values, after all, result not very different from 
those extracted from the BFKL pomeron for reasonable values of the coupling constant. 
However the infrared behavior has a strong influence on the slope value and may be essential 
for the character of the pomeron singularity. 
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Table 1. Exact and approximate gluon trajectory and po- 
tential 



q 2 


T 


T 

approx 


^>2 


-U 


approx 


0.37 


0.133 


0.138 


0.0183 


1.319 


1.293 


1.00 


0.147 


0.178 


0.0498 


1.044 


1.008 


2.72 


0.241 


0.372 


0.135 


0.695 


0.637 


7.39 


0.593 


0.767 


0.368 


0.302 


0.241 


20.1 


1.013 


1.139 


1.00 


0.0426 


0.0468 


54.6 


1.293 


1.424 


2.72 


2.78E-4 


6.77E-3 


148. 


1.515 


1.646 


7.39 


4.03E-10 


9.25E-4 


403. 


1.698 


1.829 








1097 


1.852 


1.983 








2981 


1.985 


2.116 








8103 


2.103 


2.234 









Table 1. captions 

The first column gives values of the gluon momentum squared in units A 2 . In the sec- 
ond and third columns the corresponding exact and approximate (Eq. (11)) values of the 
gluon trajectory are presented in units 6/6. The fourth column gives the values of the inter- 
gluon distance squared in units A~ 2 . In the last two columns the corresponding exact and 
approximate values of the gluon potential (with minus sign) are given (in units 3/6) 
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Table. Pomeron intercept and slope 



m 


A 


a' 


0.210 


1.02 


86.0 


0.245 


0.687 


13.8 


0.283 


0.579 


5.77 


0.4 


0.451 


1.32 


0.6 


0.372 


0.346 


0.8 


0.336 


0.150 


1.0 


0.311 


0.0812 


1.41 


0.280 


0.0333 


2.0 


0.256 


0.0138 



Table 2. captions 



The first column gives values of the gluon mass in GeV In the second column the corre- 
sponding pomeron intercepts are presented. The third column gives the pomeron slopes in 
GeV~ 2 . The QCD parameter A has been taken to be 0.2 GeV for three flavours (b = 9/4). 



